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ABSTRACT 



We analyze and model an M8.0 flare on 2005 May 13 observed by the Transition 
Region and Coronal Explorer ( TRA CE) and Reuven Ramaty High Energy Solar Spec- 
Ph ■ troscopic Imager {RHESSI) to determine the energy release rate from magnetic recon- 

nection that forms and heats numerous flare loops. The flare exhibits two ribbons in 
"^ . UV 1600 A emission. Analysis shows that the UV light curve at each flaring pixel 

I , rises impulsively within a few minutes, and decays slowly with a timescale longer than 

;h ! 10 minutes. Since the lower atmosphere (the transition region and chromosphere) re- 

t^ ' sponds to energy deposit nearly instantaneously, the rapid UV brightening is thought 

to reflect the energy release process in the newly formed flare loop rooted at the foot 
*vj ■ point. In this paper, we utilize the spatially resolved (down to 1") UV light curves and 

K^ ■ the thick-target hard X-ray emission to construct heating functions of a few thousand 

»vj . flare loops anchored at the UV foot points, and compute plasma evolution in these loops 

•/^ , using the Enthalpy-Based Thermal Evolution of Loops (EBTEL) model. The modeled 

. \ coronal temperatures and densities of these flare loops are then used to calculate coro- 

^s \ nal radiation. The computed soft X-ray spectra and light curves compare favorably 

cn ' with those observed by RHESSI and by the Geostationary Operational Environmental 

Satellite {GOES) X-ray Sensor (XRS). The time-dependent transition region differential 
emission measure (DEM) for each loop during its decay phase is also computed with 
/\ . a simplified model and used to calculate the optically-thin C iv line emission, which 

C^ , dominates the UV 1600 A bandpass during the flare. The computed C iv line emission 

decays at the same rate as observed. This study presents a method to constrain heating 
of reconnection-formed flare loops using all available observables independently, and 
provides insight into the physics of energy release and plasma heating during the flare. 
With this method, the lower limit of the total energy used to heat the flare loops in 
this event is estimated to be 1.22x10^^ ergs, of which only LOxlO^*^ ergs is carried by 
beam-driven upflows during the impulsive phase, suggesting that the coronal plasmas 
are predominantly heated in situ. 

Subject headings: Sun: activity - Sun: flares -Sun: transition region - Sun: UV radia- 
tion - Sun: X-ray radiation - magnetic reconnection 



Introduction 



Solar flares are generally believed to be a result of magnetic reconnection in the corona. After 
reconnection, the connectivity of the magnetic field changes and the field relaxes to a lower energy 
state. The energy released during reconnection is transported along reconnection-formed flare 
loops to the lower atmosphere, by conductive flux or non-thermal particles, giving rise to enhanced 
emissions in optical and ultraviolet (UV) wavelengths. As the chromosphere is impulsively heated, 
a pressure front is formed that pushes plasmas both upward (chromospheric evaporation) and 
downward (chromospheric condensation). The upflow fills post-flare loops, where energy is radiated 
away in soft X-rays (SXRs) via thermal bremsstrahlung at temperatures of up to a few tens of 
million Kelvin (MK). When these loops cool down to a few MK, they begi n to be visible in extreme 



ultraviolet (EUV). (See the comprehensive review bv lFletcher et al.ll2011 



In this qualitative picture, several fundamental questions remain unanswered in terms of the 
quantity of energy. First, we do not know how much energy is released by reconnection, and how 
much energy in total is radiated by flare plasmas (and carried away by coronal mass ejections, if 
any); in principle, the former should equal the latter. Second, we do not fully understand how, 
and by how much, flare plasmas are heated and particles are accelerated by reconnection-released 
energy, generating observed radiation signatures across the full spectral range. 

A variety of hydrodynamic simulations have been performed to understand energy transport 
in flare loops and the atmospheric respon se to energy deposition at the feet of these loop s in the 
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Since plasmas are frozen 
in magnetic loops, most of these simulations use one-dimensional (ID) models. Some of these 
hydrodynamic simulation s concentrate on the lower atmosphere by using full non-LTE forni ulation 
of ra diation transfer (i.e.. Fisher et al.lll985a ;. Hawlev fc Fisheiill994 ;_ Abbett fc Hawlevlll999 ;_ Allred 



et al. l2005l ). while others focus on the upper atmosphere doi ninated by optically thin emissio ns and 



treat the lower atmosphere as boundarv conditions (i.e.. iNagai 



1982 



1980 



Somov et al 



1981 



Peres 



Cheng et al.l 119831 : iNagai Sz Emslid Il984l ) . The radiative hydrodynamic simulations 



et al. 

have succeeded in reproducing optically thick emissions in solar flares; however, the corona could 
only reach a few MK, which is much less than observed temperatures of tens of MK. On the other 
hand, the hydrodynamic models with optically thin emissions, though unable to address emissions 
from the chromosphere, work better at reproducing high coronal temperatures and densities, and 
enhancements of UV emissions from the transition region. 

With high resolution EU V images obtained by the Transition Region and Coronal Explorer 
( TRA CE; iHandy et al.l Il999l ) , it is found that th e post-flare arcade is composed o f at least a few 
hundred flari ng loops formed succ essively (e.g., lAschwanden fc Alexander I2OOII ). It is further 
confirmed by iFletcher et al.l (|2004l ) that the flare UV ribbons are made of sma ll kernels outlining 
the feet of flare loops. To simulate sequentially formed loops in solar flares, iHori et al.l ()1997l ) 



developed a "pseudo-two-dimensional' model with multiple ID loops heated progressively from the 
innermost loop to the outermost . The multi-ID model was further used to reproduce X-ra,y and 



EUV observations of flares (e.g., iHori et al. 



Warrenll2006l : iReeves et alJl2007l : iReeves fc Moats 



1998 



Reeves &: Warrenll2002l : IWarren &: Doschek 
2O10h . 



2005 : 



The hydrodynamic evolution of plasma inside the flare loop appears to be governed by how 
the loop is heated. Specifically, we need to know when, for how long, by how much, and by what 
physical mechanism a flare loop is heated. Among existing models, many start with an unspecified 
ad hoc heating source in the corona, and the ene r gy is transported through conduct ive flux (e.g., 
Somov et al]ll982l : lNagailll98ol : Icheng et al]ll983l : iMariskalEgS?! : Isterling et a l] 1 19931) . Some other 



model s specify the source of t h e energy and/ o r mec h anism of hea. t ing. F or ex ample. ISomoy et al 



(|l98ir ). iNagai fc Emslid (119841 ). iMariska et al.1 (|l989f ). lEmslie et al.l (|l992h . and IReeves et all (120121 ) 
assumed that energy is carried by non-thermal particles and the loop is heated by Coulomb collisions 
of electrons with ambient plasmas, either in the corona (thin target) or at the lower a tm osphere 
(thick t arget), using an an alytic expression of electron energy derived by lEmslid ()1978l ) or [Brown 
(JI973I ). IReeves et al.l (J2007l ) and IReeves Sz Moatd ( 2010l) used Poynt i ng flux de rived from numerical 



solutions of the loss-of-equilibrium flare model (JLin &: Forbes 



2000 



Linll2004l ) as the energy input 



in the loop heating model, and could relate the resultant plasma radiation time profile with the 
kinematic evolution of the flux rope ejection. 



Effort has also been m ade to constrain the energy input with observations. IWarren &: Doschek 



(|2005l ) and I WarrenI ()2006l ) derived the heating rate empirically by matching the loop heating model 
results with the SXR fluxes observed by the X-ray Sensor (XRS) on the Geostationary Operational 
Environmental Satellite {GOES). In their model, the heating events are assumed t o have a trian gular 
time profile with a fixed duration (60 s in lWarren fc Doschekll2005l . and 400 s in IWarrenll2006l ) and 
their magnitude s and distributions during the flare are adjusted to match the observed SXR flux. 
Longcope et al.l (J2010l ) modeled and analyzed a flare observed by TRACE, using as the heating 
rates the reconnection-released energy, which is calculated from a patchy reconnection model. The 
model takes the observed photospheric magnetic field as the boundary, and uses the observationally 
measured reconnection rate and time distribution of reconnection-formed flare loops counted in 
EUV images obtained by TRACE. The modeled coronal radiation compares favorably with SXR 
observatio ns from GOES XRS and the Reuven Ramaty High Energy Solar Spectroscopic Imager 
{RHESSI; iLin et al.ll20o3 ). 



Along this avenue, recently, iQiu et al.l (J2012l ) have proposed that the heating rates in individual 
flare loops could be directly inferred from the time profiles of UV emission at the feet of flare loops. 
The spatially resolved UV light curves exhibit an impulsive rise, which indicates, an d therefore 



constr ains, when, for how long, and by how much a flare loop is heated. Accordingly, iQiu et al 



(|2012l ) constructed from these light curves the heating rates of 1600 loops, each with cross-sectional 
area of 1" by 1", anchored at impuls ively brightened U V pixels in a C3.2 flare observed by the 
Atmos pheric Imaging Ass embly (AIA; iLemen et al.ll2012l ) onboard the Solar Dynamic Observatory 
(5'ZJO: iPesnell et al.ll2012l ). Applying these heating rates to a zero-dimensional loop heating model. 



the E nthalpy-Based Thermal Evolution of Loops (EBTEL; iKlimchuk et al.1 120081 : ICargill et al 
2012a| ]bf) model, they computed plasma evolution in these 1600 flare loops and the synthetic coronal 



radiation in SXR and EUV passbands, which are compared with observations by GOES and AIA 
to verify the empirically determined heating rates. This method uses all available observations to 
constrain the heating rates from the input (the impulsive UV emission from the foot points) to the 
output (the X-ray and EUV emission from the coronal loops). The flare studied in their paper is 
primarily a thermal flare with little hard X-ray (HXR) emission, and ad hoc volumetric heating 
rates are used in the model. 



In this paper, we will improve the method of lQiu et al.l (|2012l ) and apply the analysis to an 
M8.0 flare on 2005 May 13 observed by TRACE, GOES, and RHESSI. The flare exhibits significant 
thick-target non-thermal HXR emission observed by RHESSI, suggestive of strong direct heating 
of the lower atmosphere that would result in chromospheric evaporation that sends energy back 
to the corona. We also include in the loop heating model this energy flux of non-thermal origin, 
using observed UV and thick-target HXR emissions as constraints. As the output of the model, we 
calculate the time-dependent 3-20 keV SXR spectrum and compare this with the observations by 
RHESSI. Finally, whereas the impulsive rise of the UV foot-point emission is considered to directly 
relate to the energy release process, the observed gradual decay of the UV emission from the same 
foot-point is governed by evolution of the overlying coronal plasma in the loop, which therefore 
provides diagnostics of the loop cooling. In this paper, we use a simplified model to compute the 
UV emission during the decay of flare loops to further compare with observations. These new steps 
help improve determination of the heating rates, and enhance our understanding of flare energy 
release, partition, and plasma heating and evolution. In the following text, we present in Section [2] 
observations of the flare, especially the UV and HXR observations, which will be used to construct 
the heating functions, and apply them to model coronal plasma evolution in Section [3l In Section [H 
we compute time-dependent SXR and UV fluxes from the model output and compare them with 
observations. Conclusions and discussions are given in Section [5l 



2. Observations and Analysis 
2.1. Overview of Observations 



In this paper, we study a GOES-class M8.0 flare that occurred on 2005 May 13 in NOAA Active 
Region 10759 located at N12E05 at the time of the flare. Figure [T] summarizes observations of the 
flare in a few wavelengths. High cadence (3 s) UV images obtained by TRAGE, given in the bottom 
panel, reveal that the M8.0 flare is a typical two ribbon flare. The flare UV ribbons expand away 
from the magnetic polarity inversion line for half an hour, indicating that magnetic reconnection 
continues to form ne w loops and relea. s e energy i n them, as depicted by the standard flare mode l 
(the CSHKP model; ICarmichaell Il964l : ISturrocklll966l : iHiravamal 1 19741 : iKopp fc: Pneumanlll976l ). 
Significant X-ray emissions are observed in this flare by RHESSI. Previous morphology studies of 



this same flare show that, during the impulsive phase, the sources of non-ther mal HXR emiss ions 
of >25 keV are located within the UV flare ribbons (see Fig. 2 and Fig. 3 of lLiu et al.l 120071 . for 
images of HXR kernels in 25-50 keV and 50-100 keV energy bands overlaid on the UV contours). 
These observations confirm that UV and HXR emissions are both produced by heating of the lower 
atmosphere (the transition region and chromosphere) by energetic particles precipitating at the feet 



of flare loops during the impuls i ve phase (see, e.g., ICheng et al. 
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The top panel in Figure [T] shows normalized and background-subtracted X-ray and UV light 
curves observed by GOES, RHESSI, and TRACE, respectively. The UV light c urve is the total 
count rate (in units of DN s~^) derived from the semi-calibrated UV images (see lQiu et al.ll2010l . 
Section 3.1 for calibration techniques), with pre-flare counts subtracted. It follows closely the 
>25 keV HXR light curve during the rise of the flare and the impulsive phase, but continues to rise 
after the HXR peak at 16:42 UT, and reaches maximum flve minutes later. Also plotted is the time 
derivative of the GOES 1-8 A SXR light curve, and it is seen that the UV light curve peaks at the 
same time (16:47 UT) as the time derivative of the SXR flux. The tem poral correlatio n between 
HXR and SXR time derivative has been known as the Neupert effect (Neupert 19681: Dennis Sz 



Zarro ll993l ). with the basic idea that non-thermal electrons precipitate at the chromosphere, losing 
their energy by Coulomb collision to give rise to HXR emission, and at the same time, driving 
chromospheric evaporation to fill post-fiare loops that are subsequently observed in SXRs. This 
fiare, however, exhibits a similar Neupert effect between the SXR derivative and UV, instead of 
HXR, light curves. This is evidence of continuous energy deposition in the lower atmosphere, most 
likely by thermal conduction, after the impulsive phase when thick-target HXR emission is no 
longer significant. 

The time sequence of UV images further shows that the continuous rise of the UV emission 
after the HXRs is produced by newly brightened UV ribbons (bottom panel in Figured]). The 
spread of the fiare UV ribbons across the longi tudinal magneti c field (lower right pan el) provides 
a measurement of magnetic reconnection fiux (JQiu et al.l l2002l : iLongcope et al.l 120071 ) , plotted in 
the top panel. The reconnection flux starts to grow at 16:30 UT at the beginning of both UV and 
X-ray light curves. The continuous increase of reconnection flux and dramatic decay of non-thermal 
emission after 16:50 UT (the post-impulsive phase) conflrm that reconnection continues to form 
new loops and release energy in them, and the thermal process is dominant in this late phase of the 
flare. The post-impulsive reconnection flux amounts to 2x10^^ Mx, about one third of the total 
reconnection flux measured for this flare. 

Comparison of the UV and X-ray observations of the flare suggests that imaging UV observa- 
tions provide information of energy release, which is then transported to the lower atmosphere by 
either non-thermal electrons or thermal conduction, in newly formed flare loops anchored at newly 
brightened UV ribbons. Based on this idea, we have analyzed spat ially resolved U V light curves, 
using them to construct heating rates of flare loops in a C3.2 flare (iQiu et al.ll2012l ). The analysis 
is applied to this flare, as well, in the following Sections. 



2.2. Characteristics of UV Light Curves 



Examples of the UV light curve in a flaring pixel (1" by 1"), shown in the top panel of 
Figure [21 typically exhibits a rapid rise, peaks within several minutes, and then decays slowly with 
a characteristic "cooling" time of more than 10 rninutes. Such characteristics have been f ound in 
UV observations of other flares as well (iQiu et alJ 120101 : ICheng et alJ 120121 : iQiu et alJl2012l ). Since 
the lower atmosphere (the transition region and chromosphere) responds to energy deposition on 
very short timescales, the rise time of UV emission reflects the timescale of energy release in the 
flare loop anchored at this pixel. The gradual decay, on the other hand, is coupled with subsequent 
cooling processes in the overlying corona. 

To characterize the rapid rise of the spatially resolved UV light curves, we flt the rise of the 
UV count rate light curve (in units of DN s~^) to a half Gaussian, 



I(t) = Jo exp 



-{t-tof 



2r2 



(t ^ to) 



fl) 



where Iq is the background-subtracted peak count rate, to is the peak time, and r is the charac- 
teristic rise time. For this flare, 5127 flaring pixels (each of size 1" by 1") are identifled from UV 
images with 3 s cadence. The lightcurve of each pixel is smoothed to 10 s and its rise phase is 
fitted to a half Gaussian. Examples of the flt to the observed light curve are plotted on the top 
panel of Figure [21 Histograms of the rise times from flts to all the flaring pixels are given in the 
bottom panel of Figure [21 showing that most flaring pixels rise over timescales from a few tens of 
seco nds to a few minutes. These rise times are svstematicallv shorter than those derived by Qiu 



et al. (|2012l ) in another C3.2 flare. This may be due to the higher cadence (3 s ) of the observations 



in this flare compared with the 30 s cadence of the other data, or more likely, the more impulsive 
nature of this flare that has significant non-thermal emission. 

Observations also show that UV pixels peaking after 16:52 UT, when the HXR emission has 
been reduced significantly, tend to evolve more slowly than those peaking earlier. This suggests 
different evolutionary timescales in both the heating and cooling phases when energy is transported 
through flare loops by different mechanisms, by non-thermal electrons or by thermal conduction. 
To examine the distinction, we compare UV light curves of pixels peaking in three stages relative 
to the HXR evolution. Pixels in the first stage (UVl) peak before 16:45 UT when HXR emission is 
significant. These pixels are marked in dark black to light blue colors in Figure [H and the sum of 
their light curves ( purple curve in Figure [2]i) has a very good temporal correlation with the HXR 
light curve. UV pixels categorized in the second stage (UV2) peak during the HXR decay from 
16:45 UT to 16:52 UT. The sum of their light curves ( orange in Figure [21i) contributes mostly to 
the peak of the total UV light curve. UV pixels in the third stage (UV3; red in Figure [211) peak 
after 16:52 UT, when >25 keV HXR emission has nearly ended. Compared with UVl and UV2, 
the UV3 light curves evidently rise more slowly and also decay more gradually. In Figure [2^, we 
plot the histograms of the rise times for pixels at the three stages, which shows that UV3 pixels 
rise more slowly on average. The differences of the UV light curves in the three stages indicates 



the different effects of beam heating and conductive heating on UV emissions, and in turn, the UV 
evolution could give us a clue to the nature of energy release, thermal or non-thermal. 

Hydrodynamic and radiative transfer models of the lower atmosphere during flares have shown 
that the lower atinosphere responds within a few seconds to the onse t of impulsive energy injection 
(JEmslie fc Nagailll985l : IPisher et alJll985al : ICanfield fc Gayleyl 119871 ). Therefore, the rapid rise of 
UV emission from the upper atmosphere or transition region may be considered to scale with the 
impulsive energy release in reconnection-formed flare loops. With this idea, we have implemented 
a method to construct heating functions and used th em to cornpute p lasma evolution in flare loops 
for a C3.2 flare, which is primarily a thermal flare (JQiu et al.ll2012l ). In this study, we apply the 
same method, using the rising UV emission at the foot-point as an indicator of energy injection 
into the newly formed flare loop (or flux tube) rooted at the foot point. We take the start-time 
of the UV brightening as the onset of the reconnection event forming the new flux tube. The rise 
time of the UV brightness gives the duration of the impulsive energy release in the newly formed 
tube, and the maximum brightness of the pixel reflects the magnitude of the energy release (or 
heating rate) in the flux tube; this simply assumes that a brighter pixel is more strongly heated. 
These observationally measured quantities may then be used to construct heating rates and study 
subsequent plasma evolution inside flare loops that are formed and heated sequentially during the 
flare. In constructing the heating functions, significant progress in the present study is inclusion of 
the heating term by non-thermal particles that are evident during the impulsive phase. 



3. Modeling Plasma Evolution in Flare Loops 

We have identified over 5000 brightened pixels of size 1" by 1" in TRACE 1600 A images, and 
assume that a half loop (or flux tube) of constant cross-section (1" by 1") is anchored at each of the 
pixels. These are half loops, because we do not identify connectivity between positive and negative 
foot-points. For each half loop, we compute the time-evolution of plasma density, temperature, and 
pressure averaged along the loop using the EBTEL model. The energy input term in the model, the 
heating rate, is constructed from observed UV count rates at the foot point. The time-dependent 
differential emission measure (DEM) is then derived from these 5000 half loops, which will be used 
to compute SXR flux and compare with observations by RHESSI and GOES in Section [H 



3.1. Loop Evolution via EBTEL 

The OD EBTEL model calculates mean properties of loop plasmas, which have been shown 
to reasonably agree with mean values from sii nulations using the ID hydrodynamic code called 
Adaptively Refined Godunov Solver (ARGOS; lAntiochos et al.Ml999l ). The OD model is highly 
efficient at computing plasma evolution for over 5000 half loops in our study of the M8.0 fiare. 
The EBTEL model solves two equations. The energy (or pressure) equation takes into account the 



prescribed heating rates {H) as the energy input term and coronal radiative loss (Re) as well as 
the total loss (Rtr) through the base of the loop (so called transition region) as energy loss terms. 
The mass equation is governed by mass flow between the transition region and the corona. In the 
EBTEL model, this flow is a result of the difference between the energy input, including the beam 
heating and conductive flux (Fq) from the corona, and the total loss (Rtr) at the base (transition 
region). During the heating phase, the energy input to the transition region dominates the loss 
term, driving upflows known as chromospheric evaporation. During the decay, the coronal plasma 
in the loop is cooling through thermal conduction and radiation; meanwhile, the loss through 
the transition region exceeds the conductive flux into it, which drives downflows, called coronal 
condensation. 

The mean coronal electron density (n) and pressure (P) in each loop evolve according to the 
EBTEL equations. 
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The mean temperature T is determined by the ideal gas law, P = 2nkBT (including both electrons 
and ions). In the equations, A;b is Boltzmann's constant, ci is the ratio of the loss through the 
transition region to the coronal radiation {Re/ Rtr), C2 is the ratio of the average coronal temper- 
ature to the apex temperature, and C3 is the ratio of the coronal base temperature to the apex 
temperature. With the symmetry assumption, EBTEL only models heating of a half loop, with 
L being the length of the half-loop. We prescribe two heating terms in the equation: the ad hoc 
volumetric heating rate, Q{t), and the energy flux carried by beam-driven upflows, T(t). The coro- 
nal radiative loss is given by Re = n^A(T), A(T) being the empirically-determ ined radiative loss 
function for optically-thin plasmas (see equation (3) of iKlimchuk et al.ll2008l. for details) . Fq is 
the conductive flux at the base of the corona, which is defined by iKlimchuk et al.l (|2008l ) as the 



location where thermal conduction changes from be ing an energy loss above to an energy source 



below. The classical form of conductive flux is used (ISpitzerlll96 
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where kq is the thermal conductivity coefficient, taken to be 1.0 x 10 ^ in cgs units. I n the EBTEL 



model , F() is saturated for large temperature gradients; in the latest version of EBTEL (ICargill et al. 
2012a| ]bf). the gravity is included in calculating of ci for semi-circular loops, while the dependence 



of C2 and C3 on gravity is negligible. 

The re are three parame ters in the EBTEL models, ci, C2, and C3. In the latest version of 
EBTEL ( Cargill et al.ll2012aH bl). which is used in the this study, ci is self-consistently determined 
by plasma properties inside the loop. Its value varies about the mean value 2.1 during flux tube 
evolution, which is not very sensitive to different heating rates in different flux tubes. Fixed values 



of C2 = 0.87 and C3 = 0.5 are used in the study. These are mean values determined from ID 
simulations. In the simulations, C2 and C3 usually change during the loop evolution, but only 
within a small range. For simplicity, in this study, we use fixed mean values of C2 and C3 for 5000 
loops throughout their evolution, considering that differences produced by using varying values of 
these parameters will become insignificant when we sum up contributions from 5000 loops. 

The critical input to the EBTEL model, also a focus of our present study, is the heating rate. 
In general, it includes two parts. The first contribution, denoted by an ad hoc volumetric heating 
rate Q (in units of erg c ni~^ s~^), is by in gifa heating in the corona. It may result from current 



dissipation, shocks (e.g., iLongcope et al 



scattering in the corona(e.g., ISomov &: Kosugi 



or even return current (e.g., iKnight &: Sturrock 



201c 



1997 



1977 



Longcope &: Guidonill201lh. electrons trapped an d 



Karlickv &: Kosugi 



Emslie 



1980 
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Holman 
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2OI2I ). In the present 



study, the exact mechanisms for this terms are not discussed. The second contribution, denoted as 
r (in units of erg cm~^ s~^), is coronal heating from the lower atmosphere due to evaporation driven 
by non-thermal electrons that precipitate at the lower atmosphere during the impulsive phase. Of 
the total flux carried by a non-thermal beam, F is a fraction of that sent back upward into the 
corona. These two terms are distinguished as they heat the corona in different ways, Q by primarily 
raising the temperature of the coronal plasma, and T by primarily raising the density - therefore the 
r term also enters the density equation. Predictably, the different heating styles result in different 
evolutionary patterns. Furthermore, these two terms play different roles in different stages of the 
flare. Specifically, the F term as driven by beams is included only during the impulsive phase when 
thick-target HXR emission is evident. 

We must note that thermal conduction alone is able to produce chromospheric evaporation, 
which, in the EBTEL model, is the consequence of coronal evolution when the conductive flux 
exceeds the loss through the lower atmosphere during the heating phase, and is therefore not 
treated as an additional coronal heating term. The F term, on the other hand, is considered in 
this study to be produced by particles that instantaneously heat the lower atmosphere, and drive 
upflow independent o/the coronal situation. Therefore, this term contributes to both energy and 
density of coronal plasmas. 

For the 2005 May 13 flare, RHESSI ohsevvations show that there is significant HXR emission 
with energy up to 300 keV during the first ten minutes of the flare. So both in situ heating in the 
corona and heating by beam driven upflow are considered in this paper. Until now it has proven 
difficult in general to distinguish these two contributions, either theoretically or observationally. In 
this study, we use an empirical method to introduce these two terms as scaled with the UV and 
HXR light curves. Determination of these two heating terms and the effect of the partition between 
the two will be discussed in the following text. 
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3.2. Constructing Heating Functions in Flare Loops 



Following iQiu et al.1 (J2012l ). we assume that the heating rate in a flare loop (or a flux tube) is 
proportional to the rise of the UV light curve at its foot-point, given as: 



Hi{t) = Qi{t)Li + Ti{t) = XI, exp 



{t - U ) 

It} 



21 



erg s cm 



; (0 < t < oo) 



(5) 



Hi is the total heating flux in the loop anchored at ith pixel, which is composed of the ad hoc heating 
rate (QiLi) and the beam-driven flux (Fj). Although we only flt the rise of the UV light curve to 
a half Gaussian, we consider the heating function to be symmetric, or a full Gaussian. Assuming 
semi-circular post-flare loops, we estimate the length of the ith half loop by Lj = f -Di, where D, is 
the distance of the foot-point to the polarity inversion line. As observed, flare ribbons expand away 
from the polarity inversion line in a rather organized manner; therefore, Lj approximately grows 
linearly with the time of flare brightening (ij). For the 5000 half loops, Lj ranges from 35 to 55 Mm. 
Although loops are in general not semi-circular, our experiments have shown that variations of the 
loop lengths within a factor of 1.5 do not signiflcantly change the synthetic total emission. 

To relate the total heating rate to UV emission, we employ a scalar A that converts the count 
rate (DN s^^) to a heating rate (erg s^^). The value of this parameter depends critically on the 
lower atmospheric response to beam heating or conductive heating and on mechanisms of UV 
emission. As a rule of thumb, we consider that A takes a larger value when conduction heating 
dominates than when beam heating dominates, for the simple reason that, with the same amount of 
energy, beam heating occurring in the lower atmosphere would generate stronger UV emission than 



conductive flux fe.g.. lEmslie &: Nagailll985l ). In this study, A takes the value of 1.9x10^ ergs DN 



for loops whose foot-point UV light curves peak before 16:48 UT, then linearly increases until 
16:52 UT; afterwards, when HXR emission has finished, A stays constant at 2.5x10^ ergs DN^^. 
Note that in this study, we do not model the lower atmospheric heating and dynamics, but instead 
use this simple empirical model to minimize the number of free parameters. The A values quoted 
here are determined by best matching the model-computed time-dependent SXR emissions with 
those observed by GOES, as will be described in the next section. These A values, combined with 
the peak count rates, correspond to the peak heating flux (H) of order 10^ to lO^'' ergs s~^ cm~^ 
for the few thousand half loops in our model. 

In the heating term, the coronal in situ heating QL is present in all flare loops; F is present only 
for loops brightened during the impulsive phase (UVl and UV2 pixels), and is gradually switched 
off after 16:52 UT (UV3 pixels). For one loop, the ad hoc heating rate and beam driven energy flux 
both have the same Gaussian time proflle of exp [— (t — ti)'^/{2T^)], and the partition of the beam 
driven energy flux, i.e., Ti/Hi, is constant. For different loops heated at different times represented 
by ti, the partition Ti/Hi is different; whereas the net heating flux Hi in a loop is proportional to 
the UV count rate at the foot-point, the beam driven flux Fj is assumed to be proportional to the 
>25 keV HXR light curve. This rather simplified treatment can be justified by the hydrodynamic 
simulations of chromospheric evaporation showing that the evaporation upflow roughly increases as 
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the heating energy flux increases ( Fisher et al.lll984l . Il985bl ) in the range of 10^-10^'' erg s ^ cm 



-2 



The estimated heating flux of the majority of flux tubes in the impulsive phase in our study is 
within this range. During the impulsive phase, the HXR spectral index does not vary significantly 
(not shown in this paper), so we can approximate the total beam energy flux, and subsequently 
the energy flux in the beam driven upflow, as proportional to the HXR count rate light curve. The 
partition Ti/Hi, therefore, is time-dependent, and is empirically given by Ti/Hi = ^mViti), where 
r/(tj) is a time-dependent function that tracks the HXR light curve, 7^ = 0.4 is the maximum 
partition of the non-thermal energy flux used in this study, and this maximum partition occurs 
when HXR emission peaks at about 16:42 UT. 

In addition to the impulsive flare heating, a constant background heating rate of order 1 x 
10~^ ergs s~^ cm~^, which is a few thousandths of the maximum heating rate constructed from 
UV observations, is imposed on each loop to produce an initial e q uilibr ium of an average coronal 



temperature of 1.8 MK and density of 6 x 10^ cm~^. iQiu et al.l ( 20121 ) have shown that neither 
the initial state of the loop nor the background heating will affect the plasma evolution as soon as 
impulsive heating occurs, since the flare heating rate is a few orders of magnitude larger than the 
background heating rate. 

The so devised heating function in a single flux tube is mostly constrained by foot-point UV 
and HXR observations, but also depends on two free parameters A and 7m. A determines the 
total amount of energy used to heat the corona, and 7^ determines the non-thermal partition. 
Since heating via upflows or via direct heating leads to different plasma evolution pattern, the 
resultant coronal radiation spectrum will differ. Our recent experiments with another flare have 
shown that model-predicted coronal radiation by high temperature plasmas is very sensitive to 
parameters defin ing the heating function but rather insensitive to other model parameters like ci 



(jQiu et al.ll2012l ). Therefore, these parameters will be eventually constrained by comparing model- 
predicted and observed X-ray emission. In this study, we start the model with an initial guess of 
the parameters, and gradually adjust them so that the synthetic SXR spectrum and light curves 
best match those observed by RHESSI and GOES. 



3.3. Evolution of Flare Plasma in One Loop 

With methods described above, we compute plasma evolution in each flux tube anchored at 
the UV flaring pixel using the EBTEL model. Figure [3] shows, in the left panel, time profiles of the 
mean temperature (solid) and density (dashed) of a single flux tube rooted at its UV foot point 
with constant cross-sectional area 1" by 1". It is heated impulsively by a Gaussian-profile heating 
rate (red dot-dashed line in the right panel) constructed from the UV light curve (black solid line 
in the right panel) with r = 50 s, and Hmax = 9 x 10^ ergs s~^ cm~^, of which 40% is carried by 
the beam-driven upflow. The A and 7^ values used here are determined from comparing synthetic 
SXR light curves and spectra with observations, which will be described in the following sections. 
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Coronal temperature rises as the flux tube in the corona is heated by ad hoc heating. At 
the same time, the lower atmosphere is heated by non-thermal electron beams and conductive 
flux, giving rise to enhanced UV emission as observed by TRACE. Chromospheric evaporation is 
driven to fill the coronal loop. An increase in coronal density leads to enhancement of coronal 
radiation (red dashed line in the right panel), which then cools the plasma in the flux tube. The 
coronal temperature begins to decline immediately after the peak of the heating, but the density 
continues to grow until the energy loss in the transition region exceeds the conduction flux, causing 
a downward flux, or coronal condensation. The coronal density and radiation then decline, and the 
flux tube experiences a long decay. Notably, the observed UV 1600 A emission decays on the same 
timescale of coronal evolution. This has been previously known as optically-thin transition region 
lines, such as the resonant C IV line that dominates the 1600 A b roadband UV emission, behave as 



a "coronal pressure gauge" (|Fisheiill987l : iHawlev fc Fisheiill992l ). 



3.4. Effect of Beam Driven Flux on Plasma Evolution 

To determine 7m, we investigate the effect of beam driven flux in the EBTEL model. We 
examine the evolution of coronal plasma in the same flux tube shown in Figure [3] with the same 
impulsive heating rate Hmax = 9 x 10^ ergs s~^ cm~^ but with varying partition, i.e., V /H values. 
The time profiles of coronal- averaged temperature, density, and pressure with different T /H values 
are plotted in Figure HI 

It is shown that, in general, a greater amount of beam driven flux results in lower temperature 
and higher density. If the beam driven flux does not dominate, i.e., V /H < 0.5, the average 
temperature decays to 10 MK in a few minutes, and then continues to decay toward the pre-flare 
value. From about ten minutes after the heating, the decay of coronal temperature and density 
is nearly identical for different T/H values. The same evolutio n pattern of the apex densit y and 



temperature is also displayed in ID hydrodynamic simulations by lWinebarger fc: WarrenI (|2004l ). who 
modeled loop heating with thermal conduction. Their study showed that a flare loop impulsively 
heated by the same amount of the total energy, but with varying magnitude and duration and 
at different locations along the loop, would reach the same equilibrium point when radiation and 
conduction are comparable, and the evolution of the apex density and temperature is identical 
thereafter. 

On the other hand, when the beam driven flux dominates the energy budget with T/H > 0.5, 
the flux tube will attain much higher density in a short time due to strong upflow. Meanwhile, 
less direct heating Q results in a lower peak temperature. The two effects, high density and low 
temperature, lead to faster radiative cooling in the decay phase. Therefore, the flux tube evolves on 
much shorter timescales. A similar evolution pattern has also been produce d in the ID numeri cal 



simulation that takes into account beam heating with different pitch angles (JReeves et al.ll2012l ). 



The experiment above suggests that varying the partition parameter T/H will lead to different 
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temperature and density in the heating phase, and for very large T /H values, evolution of the flux 
tube in the decay phase is also modified significantly. Therefore, the experiment provides us with 
two ways to estimate the optimal T /H value in this flare. First, we may compare the synthetic SXR 
spectrum using the DEM from modeling over 5000 flux tubes with the X-ray spectrum observed 
by RHESSI during the heating phase. Figure [5] shows the DEM constructed using different 7^ 
values for the 5000 tubes. It is seen that, given the same amount of total energy, the larger value 
of 7m, i.e., the larger amount of energy carried by beam driven upflow, will lead to greater DEM 
increase at lower temperatures, whereas the DEM at the high temperature end would be reduced. 
Such variation will affect the synthetic SXR spectrum formed at the temperature range of a few to 
a few tens of MK. Therefore, the value of 7^ can be estimated by comparing the synthetic SXR 
spectra with RHESSI observations, as will be presented in Section 14.21 



Second, we will compare the observed and modeled timescales of the foot-point UV emission 
during the decay phase of the flare, based o n the principle that UV lirie emi s sion behaves as a 



coronal pressure gauge in the decay phase (JFisher et al.l Il985al : iFished 119871 : iHawley &: Fisher 



I992I ). These will be elaborated on in Section 14.31 As a quick look, the right panel of Figured] 
shows the observed UV light curve versus pressure; it appears that the UV emission decays on the 
same timescale as the coronal pressure for T /H < 0.5. 



3.5. Properties of Coronal Plasmas in Multiple Flare Loops 

Our methods yield the best estimate of A and 7m, with values given in Section [3.21 with which 
we can determine heating functions of over 5000 flux tubes. The left panel of Figure [6] shows the 
time profiles of the total in situ heating rate (^ QiLi) and total beam heating rate (^ Fj) summed 
for all the flux tubes, in comparison with the observed UV total counts light curve and HXR 25- 
50 keV counts light curve. The total heating rate (^ Hi) is shown nearly scaled with the total 
UV count rate light curve from the rise to the peak - note that the decay of the UV light curve 
is governed by coronal evolution, which is not part of the heating. The beam driven energy flux is 
proportional to the observed >25 keV HXR light curve. The time integral of the total heating rate 
and beam heating rate yields the estimate of the total energy used to heat the corona and the total 
energy carried by beam-driven upflows in this flare, which are 1.22 x 10^^ ergs and 1.9 x 10^'^ ergs, 
respectively. The figure shows that the largest non-thermal partition amounts to about 40%, which 
occurs around the peak of HXR emission at 16:42 UT. Note that although the total heating flux 
Hi in a single flux tube is considered to be proportional to the UV light curve at the foot-point, 
the sum of the total heating fluxes in all flux tubes, ^ ff « = Y^{QiLi + Ti), is scaled with the total 
UV light curve by different values of A during different phases of the flare, as can be seen from 
the Figure. A is time-dependent in our method; given the same a mount of energy, beam heating 



produces greater UV emission compared with conduction heating (JEmslie &: Nagailll985l ). 



With these heating rates, we compute plasma evolution in over 5000 flare loops formed by 
reconnection and heated at different times and by different amounts of energy. The distribution of 
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the peak temperatures and densities of these flux tubes is shown in the right panels of Figure [6l 
For this flare, the peak temperature of most flux tubes varies from 7 to 26 MK, and the peak 
temperature distribution appears bi-modal. The second peak in the distribution at around 22 MK is 
mostly contributed by flare loops formed around 16:45 UT, and is close to the effective temperature 
derived from the ratio of th e two-channel (1-8 A and 0.5-4 A) emission measured by GOES with 
the isothermal assumption (jWhite et alJl2005l ). The peak density of the flux tubes ranges from 
5-30x10^ cm-3. 

The DEM of the coronal plasma is readily calculated from the temperatures and densities of 
these loops. Figure [7] shows the time evolution of the coronal plasma DEM averaged every 20 
seconds for the M8.0 flare. Note that the higher DEM at about 2 MK at the beginning of the 
flare is due to background heating. It is seen that in the rise phase of the flare, the DEM increases 
toward higher temperature, suggesting that more flux tubes are heated to higher temperatures, 
while during the decay phase, the DEM rises at lower temperature, reflecting cooling of the flux 
tubes. The time evolution of the coronal DEM will be used to calculate time-dependent SXR 
emission to verify the model results, as discussed in the next Section. 



4. Comparison with Observations 

The time-dependent coronal DEM computed from over 5000 flux tubes is convolved with instru- 
ment response functions to compute synthetic SXR light curves and spectra, which are compared 
with those observed by GOES and RHESSI (Sections 14.11 and 14.21 respectiv ely). The transition 
region DEM during the decay phase is also derived from the coronal pressure (JFisher et al.lll985al : 



Fisheiill987l : iHawley &: Fisherlll992r). with which we calculate th e the C iv line emission using the 
CHIANTI atomic database (jDere et al.1 119971 : Landi et al.ll2012l ). and compare it with UV 1600 A 
observations (Section! 



4.1. Comparison with GOES Light Curves 

The two X-ray Sensor (XRS) photometers onboard GOES measure full-disk integrated SXR 
emissions in two energy bands, 1-8 A and .5-4 A, with 3-second c adence. The solar 1-8 A flux, 
which is dominated by continuum emission (|Mewdll972l : lKatdll976l ). is believed to originate from 
hot (T > 10^ K) plasma in c oronal loops in active regions. T he full-disk SXR images observed by 



the Solar X-ray Imager (SXI; iHill et al. 



2005 



Pizzo et al.ll2005l ) onboard GOES indicate that NOAA 



Active Region 10759, where the M8.0 flare occurred, is the dominant SXR source on 2005 May 
13 from 16:00 UT to 20:00 UT (images not shown here). Therefore, the background-subtracted 
GOES light curves reflect SXR emissions by the M8.0 flare. The comparison of the synthetic SXR 
light curves from the modeled flare plasma with the observations by GOES XRS yields the best-flt 
model parameter A. 
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Figure [8] shows the synthetic SXR emissions in two channels plotted against the observed light 
curves. In both band-passes, the synthetic light curves follow the observed ones very well from the 
rise until 16:50 UT. In the 0.5-4 A channel, the modeled and observed fluxes start to decline at 
around the same time, but the modeled flux decays more rapidly than the observed. In the 1-8 A 
channel, the modeled radiation flux begins to drop while the observed radiation continues to rise for 
another 5 minutes. We also note that at the start of the flare before 16:40 UT, the observed GOES 



1-8 A flux starts to rise earlier than the modeled flux. As discussed by lQiu et al.l ( 20121 ). simply 
changing model parameters cannot compensate for the flux deficiency; instead, these discrepancies 
are most likely caused by weak heating events in the very early and late phases of the flare, which 
might not be identified in the UV foot-point emission. 



4.2. Comparison with RHESSI Soft X-ray Spectra and Light Curves 

RHESSI is designed to observe solar high-energy emissions from SXRs to gamma-rays (3 keV 
up to 17 MeV) with an unprecedented combination of high time, spatial, and spectral resolutions 
(JLin et al.ll2002l ). Whereas the thick-target HXR (>20 keV) observations provide us with a guide 
for constructing heating functions, we will further compare the model computed SXR (<20 keV) 
spectrum with that observed by RHESSI. The 3 to <20 keV X-ray emission is usually thermal 
bremsstrahlung radiation produced by plasmas in the flare loops heated to temperatures of a few 
to a few tens of MK. RHESSI images confirm that this emission is generated in the corona. 

We calculate the time-dependent SXR spectrum using the DEM distribution obtained from our 
model (see Figure [7j). At each time (20 s cadence), we first calculate the X-ray spectrum observed 
at Earth (in units of phot ons cm~^ s~^ keV"^) from the opt ically-thin thermal bremsstrahlung 



radiation using CHIANTI (JDere et al.l 119971 : iLandi et al.l l2012l^. including both th e line and con- 
tinuum cont ributions and assumin g solar coronal abundances (JFeldman et al.lll992l ) and ionization 



equilibrium (iMazzotta et al.lll998f). The calculated spectrum is then convolved with the RHESSI 



detector response (jSmith et al.l 120021 '). obtained from pre- flight and in-flight instrument calibration 
and modeling, to convert the photon spectrum to the count spectrum that would be directly ob- 
served by RHESSI. The detector response accounts for instrumental effects due to pulse pileup, 
escape of K-shell fluorescence photons from the detectors, attenuation by the aluminum shutters, 
and the ene r gy-de pendent detector effective area, and includes the additional in-flight corrections 
from ICaspil (J2010l ). To get the best comparison of the sp ectra, only detecto r G4 is used since it 
had the best in-flight resolution at the time of the flare (jSmith et al.ll2002l ). By integrating the 
synthetic spectrum over different photon energy ranges, we also obtain synthetic SXR light curves 
in units of counts rate. 

Figure [9] shows the synthetic SXR light curves in the 3-6, 6-12, and 12-25 keV bands, respec- 
tively, in comparison with RHESSI observations. For most of the impulsive phase, the synthetic 
light curves agree very well with the observed ones. The synthetic fluxes are lower than observed at 
the start of the flare, for a few minutes after 16:37 UT, when RHESSI emerges from eclipse. They 



16 



then rise rapidly to catch up with the observed fluxes as more flare loops formed and heated by 
reconnection are identified in the UV foot-point observation. The computed and observed fluxes 
overlap for 15 minutes from 16:44 UT until 16:59 UT, well after the impulsive phase. Afterwards, the 
computed fluxes start to drop faster than observed. This is consistent with the model-observation 
comparison of the GOES li ght curves, which u s ually yield slightly lower plasma temperatures than 



RHESSI light curves fe.g.. iHannah et al.ll2nn8l : ICaspi fc Linll2nim . 



More details can be learned by comparing the SXR spectra. Figure [TOl illustrates the synthetic 
and observed spectra at a few different times, indicated in the top panel of Figure [U during the 
flare. (For a complete view, a movie comparing the modeled and observed SXR spectra, from 
16:37 to 17:15 UT with 20 s cadence, is available online.) Figure [TT] summarizes the comparison of 
the time-dependent spectra with two quantities, the ratio of the synthetic to the observed counts, 
p{e) = Cm{()/Co{e), averaged over the photon energy range e from 6 to 15 keV, and the slope 
a of this energy dependent ratio versus the photon energy, which is obtained by fitting the ratio 
to a linear function of the photon energy p{e) = ae + h. If the synthetic spectrum is the same 
as the observations, p = b = \ and a = 0. p is a measure of magnitude comparison, and we 
regard empirically that the synthetic and observed spectra agree with each other when p ranges 
between 0.7 and 1.3. a is an indication of the plasma temperature distribution, since the thermal 
bremsstrahlung spectrum is temperature dependent. When a > 0, the model has produced more 
high temperature plasma than observed, and when a < 0, there is a lack of hot plasma in the 
modeled DEM. 

In these plots, we only compare the spectrum up to 20 keV, beyond which the spectrum is likely 
dominated by thick-target bremsstrahlung produced by non-thermal electrons colliding at the lower 
atmosphere. For most of the impulsive phase, from 16:44 UT to 16:59 UT, the synthetic spectra 
agree with the observations (see panels b, c, and d in Figure [T0|) . with mean p value between 0.7 
and 1.3 and a close to zero. From 16:50 to 16:55 UT, around the maxima of the GOES SXR light 
curves, the modeled thermal emission dominates up to 20 keV (see panels c and d), indicating strong 
heating in flare loops although the non-thermal HXR emission has become insignificant. These plots 
showing good agreement between the model and observations suggest that the model has rather 
accurately reproduced the DEM responsible for the thermal bremsstrahlung X-ray radiation in the 
observed range during this period. 

From the start of the flare until 16:42 UT, the synthetic spectrum is lower than observed by 
nearly the same fraction across the energy range up to 15 keV (panel a in Figure [T0|) . The synthetic 
flux at >15 keV is still smaller than observed, which may indicate a lack of high-temperature plasma 
in the model, but may also be due to predominant non-thermal bremsstrahlung radiation down 
to 15 keV. RHESSI images in the 15-20 keV energy range during this period (not shown) show 
that the X-ray source is rather extended, likely including emissions from both the foot-point and 
loop-top, which cannot be resolved with good photometric accuracy. Therefore, both scenarios are 
plausible. In the late phase of the flare after 17:00 UT, the modeled flux drops below the observed 
flux again. In this phase, thick-target HXR emission has ended, so the low flux above 15 keV is 
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likely caused by a lack of high temperature plasma in the model. We also note that in the early 
and late phases of the flare, when the average ratio p is far below 1, the slope a becomes negative. 
These indicate that the model does not generate enough high temperature plasma during these 
periods, maybe due to a lack of heating events identified in our method. The bottom panel of the 
figure shows the number of heating events identified during the flare. It is seen that the number 
of heating events drops rather significantly both before 16:39 UT and after 16:55 UT. It appears 
that the shortage of heating events leads to insufficient X-ray flux in the calculation about 3 to 5 
minutes later, which is roughly the timescale for temperature and density increases in newly formed 
flux tubes. The lack of heating events identified from lower atmospheric signatures (brightened UV 
or HXR foot-points) could be explained by the fact that there is neither strong conduction flux 
nor significant non-thermal beam to deposit enough energy at the lower atmosphere to generate 
observable signatures. In the rise phase, it is plausible despite the significant non-thermal HXR 
emission, if the non-thermal electrons are partially or largely trapped in the corona and do not 
reach the chromosphere to generate HXR foot points; the trapped electrons would then be a source 
of in situ heating in the coron a. Such a mechanisin has been proposed for the "pre-impulsive" 
phase of a few large flares (e.g., iLin et al.ll2003l : ICaspill201Cl : ICaspi fc Linll201Ci ) where non-thermal 
emission is observed in the corona with no identifiable HXR foot points, and may also contribute 
here. 



Comparison between the observed and synthetic spectra, as characterized by p and a, provides 
us with observational constraints to the heating function parameters A and 7^, in addition to the 
observed UV and HXR signatures at the flare foot-points that constrain the time profiles of the 
heating rate. 



4.3. UV Emission in the Decay Phase 

The long decav of TRACE 1600 A emission in individual flaring pixels is reported bv Qiu 



et al. ( 2OI0I ) and lCheng et al.l ( 20121 ). Flare emission in this broadband channel includes both the 
enhanced UV continuum and the C iv line emission. The C IV resonance doublet at 1548 A and 
1551 A could only be emitted over a very narrow temperature range around 10^ K, and therefore the 
C IV intensities rise when the transition region is heated with more plasma raised to 10^ K during 
the flare. The C iv irradiance was observed to be enhanced by three to four orders of ra agnitude 
over the pre-flare emission during the impulsive phase of a large flare (JBrekke et al.lll996r). During 
the decay phase, the C iv emission is domi nated by the evolution of coronal pressure ()Fisherlll987 : 
Hawlev fc FisheJl992l : [Griffiths et al.lll998h . The solar UV continuum radiation in the 1000-2000 A 



wave l ength range is also enhanced during the flare (e.g., ICook h BruecknerJll979l : ICheng et al. 
19841 . Il988l ). In particular, the continuum emission below 1682 A is prim arily contributed by bound- 
free transitions of Si I, w hich is primarily excited by the C iv doublet (JMachado &: Henouxlll982 : 
Machado &: Mauaall986l ). Therefore, th e continuum intensit y at A < 1682 A is approximately 
proportional to the C IV line intensity (jPhillips et al.lll992l ). Though there is no direct C iv 
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measurement for the M8.0 flare studied in this paper, observations of stehar flares show that the 
time profile of C iv is similar to the U V 1600 A emissioii in solar flares observed by TRACE, both 
exhibiting a fast rise and a long decay (jVilhu et al.lll998l ). 



In this study, we have used the rapid rise of UV emission to constrain the time, duration, and 
magnitude of the heating rate in individual flare loops. The gradual decay of UV emission from 
the same foot-points reflect the evolution of the over-lying corona that has been heated impulsively 
and then cools down on a much longer timescale. Using the coronal plasma properties from the 
model, we further calculate the C iv line emission during the decay phase to compare with the 
TRACE UV observation. 

To calculate the C iv line emission o i i the solar surface, it is important to know the opacity. 
Doschek et al.l (|l99ll ) and iDere &: MasonI (|l993l ) measured the intensity ratio of two C iv lines at 
A1548 and A1550 ob served on the d isk, and found the ratio very close to 2:1, which is expected for 
optically thin lines. iDoschekl (|l997l ) also estimated the opacity of C iv to be 0.099 in the active- 
region spectrum. We therefore consider that the C iv line is optically thin. The C iv photon flux 
at the surface is then calculated through 



/c jy = Ab / C(r, n) DEM(r) dT photons cm" 



s sr 



(6) 



where ^4^ is the abundance of Carbon relative to Hydrogen. In this pape r the solar coronal abun - 
dances of iFeldman et al.l (jl992l ) are used; if photospheric abundances (i.e.. lGrevesse fc Sauva]|ll998l ) 
are used, the C IV line emission will be red uced by 15%. C(T, n) is the contribu tion function for the 



C IV hue and is obtained from CHIANTI (Dere et al. 



1997 



Landi et al.ll2012l ) with the ionization 



equilibrium of iMazzotta et al.l (119981 ). It does not change significantly when the electron density n 
varies from 10^ cm~^ to 10^^ cm~^. DEM(r) = n? (^) is the transition region DEM, where s 
measures distance along the flux tube from the base of the transition region. For each flux tube, 
we assume that the transition region is nearly in static equilibrium in the decay phase when the 
heating has finished, which is roughly when the coronal density peaks. In this case, thermal con- 
duction is balanced by radiation, and pressure is uniform from the corona through t he transition 
region when gravity is ignored, and the transition region DEM can be calculated as (lFisheiill987 : 
Hawlev fc Fisheilll992l ): 



DEM(r) = P 



(7) 



where C(^) = fx T'^''^A{T') dT', Tq is the temperature at the base of the transition region, taken 
to be a nominal 10^ K, and A(T) is the optically-thin radiative loss function. 

The C IV line emission of one flux tube in its decay phase is calculated with equations ^ and 
([Tj) and convolved with the TRACE 1600 A band response function to get synthetic UV flux in 
units of DN s^^. The left panel of Figure [T2] shows an example of the computed C iv line emission 
in one flux tube. The computed C iv flux appears to evolve along with the observed decay, though 
the synthetic flux is smaller than the observed by a factor of 2. The right panel of the figure shows 
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the sum of the synthetic UV flux in all flaring pixels in comparison with the observed total UV 
emission in the decay phase. The observed UV decay in a pixel is obtained by subtracting from 
the UV light curve the full Gaussian derived by fitting the UV rise so that the transition region 
response to the heating is excluded. Similar to the single-pixel comparison, the computed total C iv 
emission decays on the same timescale as observed, and the computed flux is lower than observed 
by a factor of 2 to 3. 

We find this comparison satisfactory for the following reasons. First, the transition region 
is not exactly in static equilibrium. During the decay phase, the transition region loss usually 
exceeds the conductive flux, which leads to downflows from the corona, and therefore the coronal 
density gradually decreases. The downflow is literally an energy input into the transition region 
and gives rise to a greater transition region DEM than in the case of static equilibrium. Second, 
the UV continuum in the decay phase may be also enhanced due to irradiation from the C iv 
line. The observed UV emission is the sum of both the line and continuum emissions, though 
the continuum intensity is considered to be proportional to the line intensity. All these reasons 
would explain the observed UV emission being larger than the computed C iv line emission with 
the static equilibrium assumption. It is, nevertheless, very encouraging to flnd that the computed 
line emission decays at the same rate as the observed UV emission. In contrast, the computed 
X-ray emission decays faster than observed, as shown in the previous Section. Since the UV decay 
is derived from the same pixels whose rise phase is used to construct the loop heating rate, such 
comparison indirectly conflrms that the rapid decay of the computed SXR flux most likely indicates 
the presence of additional heating events in the decay phase, which are not found in foot-point 
UV radiation signatures. This is consistent with decades of observations which have shown that 
coronal plasma temperatures decrease mo re slowly than wou l d be expected from s imple calculations 
of radiative and conductive cooling (e.g., iMoore et al.lll980l : IVeronig et al.ll2002al jbl). implying that 



additional heating must be taking place despite an apparent lack of accelerated particles. These 
additional heating events may thus result from continuous reconnection in the corona that provides 
direct heating without signiflcantly accelerating non-thermal particles, or which accelerates particles 
sufficiently weakly such that they thermalize in the ambient plasma at the loop-top; both scenarios 
lead to increased SXR coronal emission without associated chromospheric emission that could be 
identified in Ti?^ C-E broadband UV or RHESSIUXR observations. A similar mechanism has been 
propo sed for the so-called "EUV late phase" emission from eruptive (and often CME-associated) 
fiares ( Woods et al.ll201ll ). where reconnection high in the corona heats post-fiare loops, resulting 



in increased EUV coronal emission without associated foot-point brightenings. 



5. Discussion and Conclusions 

We have analyzed and modeled an M8.0 flare on 2005 May 13 observed by TRACE, RHESSI, 
and GOES, in order to determine heating rates in a few thousand flare loops formed by reconnection 
and subsequently heated. The rapid rise of spatially resolved UV emission from the lower atmo- 



20 



sphere is considered to be the instantaneous response to heating in the flare loops rooted at the foot 
points, and therefore is used to construct heating rate s in individual fl are loops to describe when, 
for how long, and by how much flare loops are heated ( Qiu et al.ll2012l ). With these heating rates. 



we compute the coronal plasma evolution in these loops using the OD EBTEL model, and compare, 
for the first time, the computed time-dependent coronal radiation with observations by RHESSI 
and GOES, and as well the computed UV line emission with Ti?^ Ci? observations during the decay 
phase of flare loops. The comparison constrains parameters of the loop heating rates, which define 
the amount of total energy fiux and the fraction of energy flux carried by chromospheric upflows 
driven by non-thermal beams. 

The M8.0 flare studied in this paper is observed to have signiflcant thick-target HXR emission, 
suggesting that the lower atmosphere is also heated by electron beams during the impulsive phase. 
In this process, chromospheric evaporation is driven carrying mass and energy flux back to the 
corona. In this study, we experimentally scale the beam-driven energy flux with the observed HXR 
light curve, and examine how varying the fraction of beam heating changes the modeled coronal 
temperature and density evolution. In general, for a given total amount of heating flux, a larger 
fraction of energy carried by beam driven upflow leads to higher density, lower temperature, and 
more rapid evolution of the coronal plasma. These properties affect the computed SXR spectrum 
and its time evolution, and can be compared against the spectrum observed by RHESSI. Our 
experiment yields an optimal set of parameters, with which, the model computed SXR light curves 
and spectra agree very well with observations during the rise phase and peak of the flare for more 
than 10 minutes. 

Apart from comparison with the SXR observations, we also compute the C iv line emission at 
the foot of the flare loop during the decay, whose rise phase is used to construct the heating rate in 
the same loop. Emission in this optically-thin line during the decay phase is governed by plasma 
evolution in the overlying coronal loop and contributes signiflcantly to the UV emission observed 
in the TRACE 1600 A bandpass. It is shown that the computed C iv flux is about one third of the 
observed UV flux, which also includes the continuum radiation, and decays on the same timescale 
as observed. This experiment shows an avenue to model heating and cooling of spatially resolved 
flare loops with self-consistent constraints by high-resolution UV observations from the input (UV 
rise) to the output (UV decay) in the same flare loop. 

These results suggest that our method, which employs all available observations to constrain 
flare loop heating model, is able to capture the distribution of impulsive heating rates in numer- 
ous flare loops formed and heated sequentially by magnetic reconnection. The optimal heating 
function parameters determined for this flare yield the peak heating flux in the range of 10^- 
]^q10 gj.gg gjji-2 g-i £qj. Qygj. 5000 flare loops of cross-sectional area 1" by 1". Notably, only <40% 
of this energy is carried by beam-driven upflows during the impulsive phase, with the remaining 
>60% contributed by in situ heating in the corona. The beam heating occurs mostly during the 
first ten minutes, when the flare exhibits significant thick-target HXR emission; meanwhile, the in 
situ heating energy in the impulsive phase amounts to 6.2x10^'^ ergs, over 3 times greater than the 
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beam-driven contribution. In the post-impulsive phase (after 16:50 UT), the continuously expand- 
ing UV ribbons and the observed Neupert effect between the SXR derivative and UV (rather than 
HXR) hght curves indicate that new flare loops are formed by continuous reconnection, and in situ 
heating in the corona is predominant, amounting to 4.1x10^'^ ergs out of the total, with negligible 
beam heating. The total energy used in coronal heating in this flare amounts to 1.22x10^^ ergs, 
of which the total energy carried by beam-driven upflows is estimated to be only 1.9x10'^'^ ergs 
(~16%). Therefore, it appears that in situ coronal heating and thermal conduction play an im- 
portant role in the energetics of this M8.0 flare, in contrast to commonly accepted models that 
consider the coronal thermal plasma to be largely a by-product of beam-driv en heating of the lower 



atmosphere (chromosphere and transition reg ion) (see also discussions by iLongcope et al. 



Caspi fc Linll2010l : iLongcope fc Guidonill201ll ). A survey of 37 M- and X-class flares (jCaspi et al. 



201C 



20131 ) suggests that in situ heating may be signiflcant even down to mid-C class, and therefore this 



heating mechanism must be considered in future studies of flare energetics. 

We note that these energies set lower limits for this M8.0 flare. First, despite the excellent 
agreement between the model and observation in the rise phase of the flare, the model computed 
X-ray flux decays more quickly than observed by both GOES and RHESSI. It is most likely that 
magnetic reconnection and energy release continue in the high corona during the decay of the flare. 
These heating events might not produce s ignificant UV foot-point emissi on, and are therefore not 



identified in our method. For comparison, iKazachenko et al.l (J2009l . |2012| ) analyzed the same flare 



and estimated the total energy from GOES (including energy loss from radiation, conductive cooling 
and enthalpy flux) to be 3.1x10^^ ergs. Second, in our experiment, we do not model heating and 
dynamics of the lower atmosphere, but use a simple scaling relationship to estimate the energy flux 
carried by non-thermal beam-driven upflows back to the corona. This amount of energy is a fraction 
of the total energy carried by non-thermal beams; the remainder is lost in the lower atmosphere. 
As a reference, we fit the HXR spectrum to a power-law distribution with the standard RHESSI 
software package (not shown in the paper), and, with the low-energy cutoff and spectral index from 
the fit, estimate the total non-thermal energy to be 7.6x10'^'^ ergs, about 4 times the energy carried 
in the upfiow but, interestingly, still somewhat smaller than the energy required to heat the coronal 
plasma ■- though we note that the non-thermal energy estimate is only a loose lower bound since it 
depends critically on the fit low-energy cutoff, which is unbounded from below due to obscuration 
by the highly-dominant thermal emission. 

Our experiment provides a novel method for investigating energy release in solar flares, which 
is governed by reconnection and substantiated in formation and heating of numerous flare loops. 
This approach has a few advantages. First, by analyzing the foot-point UV signatures combined 
with HXR observations, for the first time, we are able to identify and characterize the energy release 
process in a few thousand flare loops down to 1" by 1" scale, which is thought to be close to the basic 
scale of flux tubes formed in patchy reconnection. Using the same UV data, we also measure the rate 
of magnetic reconnection and hence are able to establish the relationship between reconnection and 
energy release in temporally and spatially resolved manner. Second, using the OD EBTEL model. 
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we are able to efficiently compute plasma properties of this large number of flux tubes, and naturally 
generate a time-dependent differential emission measure (DEM) from these flare loops formed and 
heated at different times and evolving independent of one another. Subsequently, we can predict 
X-ray and UV radiation signatures directly comparable with observations, which provides further 
constraints to our determination of energy release rates. 

The method utilizes emission signatures at all available wavelengths in both the foot-points and 
overlying coronal loops. The dynamics of the lower atmosphere and corona are strongly coupled; 
they are governed by different yet coherent physics, and hence should be studied coherently. In 
this study, based on the insight from previous theoretical and numerical research, we employ some 
empirical laws to prescribe the relation between heating rates and radiation signatures. The optimal 
parameters obtained from the analysis can provide a reference for further investigation based on 
physical models, such as models of lower atmosphere heating and dynamics. We also recognize that 
the OD EBTEL model has limitations in accurately describing coronal plasma properties along the 
loop. Nevertheless, the model is highly efficient in dealing with a few thousand flare loops and 
yields the first-o rder heat i ng ra t es as useful inpu t s for mor e sophisticated loop heating models 
such as those by iMariskal (jl987l ) , lAntiochos et al.l (|l999l ) , or I Winter et al.l (J201ll ) . Our study of 
this event shows very good agreement between observed and modeled high temperature plasma 
radiations (e.g., by RHESSI and GOES), but lacks observations of and therefore comparison with 
low temperature plasma signatures, which are crucial in understanding the flare loop evolution 
beyond the conduction-dominant regime. The recently launched SDO has observed many flares 



with AIA, the high-resolution imaging telescope, in a number of UV and EUV bands (JLemen et al 



2013), and as well with the br oad-band, high reso lution spectrometers on the Extreme Ultraviolet 
Variability Experiment (EVE: IWoods et al.ll2012l ). These new observations, which cover coronal 
temperatures from <1 MK up to ~20 MK, can be combined with RHESSI observations that provide 
information of hot (>10~50 MK) plasmas, as well as information about the non-thermal emission, to 
fully observationally characterize the temperature distribution and its evolution (Caspi et al. 2013b, 
in prep.), which can be compared with - and used to constrain - our modeling. These new facilities 
will allow us to analyze and model flare plasmas following the entire process of heating and cooling, 
and also to study resolved individual flare loops in multiple wavelengths in more details, which will 
improve determination of the energy release rate during magnetic reconnection. 

Finally, the poor comparison between model and observations during the decay phase of the 
flare reflects some limitations of the present method. These can be addressed by more sophisticated 
one-dimensional modeling that also includes the lower atmosphere, to help determine the amount 
of lower atmosphere radiation, such as the UV 1600 A emission used in this study, as dependent 
on the heating mechanism and the amount of heating. Such effort will help clarify whether the 
shortage of computed flux is produced by the method missing weak heating events, or by other 
significant effects related to imperfect model assumptions. 
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Fig. 1. — Light curves and images of the M8.0 flare observed on 2005 May 13. The upper 
panel shows background-subtracted and normalized light curves in HXRs (25-50 keV) observed 
by RHESSI detector 4 (blue), SXRs (1-8 A) by GOES XRS (red), and UV (1600 A band) by 
TRACE (dark green). Also plotted are the time derivatives of the 1-8 A SXR flux (brown), and 
magnetic reconnection flux (black) measured by using UV 1600 A observations and the longitudinal 
magnetogram from SOHO/MI}!. The left and middle images in the lower panel are snapshots of 
UV images during the impulsive and decay phases of the flare observed by TRA CE, while the lower 
right panel shows evolution of the UV brightening on top of the MDI longitudinal magnetogram 
taken at 16:03:02 UT with the color bar indicating the start times of UV brightening at different 
locations. 
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Fig. 2. — Upper: observed UV 1600 A light curves (thick black line) of three different flaring 
pixels, superimposed with the Gaussian function (thin red line) which fits the impulsive rise of each 
individual light curve. Bottom left: normalized background-subtracted UV light curve of all flaring 
pixels, and UV light curve of pixels peaking at three different stages (see text), compared with the 
HXR 25-50 keV light curve. The UV light curves are normalized to the maximum of the total UV 
light curve. Bottom right: histograms of the rise times of UV light curves for pixels brightened at 
three different stages as shown in the bottom left panel. 
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Fig. 3. — Evolution of plasma properties in a flare loop with cross section 1" by 1" rooted at one 
flaring pixel. Left: computed time profiles of coronal-averaged temperature (solid) and density 
(dashed) of the flare loop. Right: observed UV 1600 A count rate light curve (solid black line), the 
constructed heating function (dot-dashed red line), and computed coronal radiation rate (dashed 
red line) of the loop. 



25 


(a) H 


r/H - 


1 20 






0.1 . 

0.2 


i 15 




\ 


0.3 

0.4 ■ 


P 




\\ 


0.5 


Q- 10 




0.6 : 

0.7 


5 




\^ 


^^"- 


n 







20 



30 40 50 60 70 80 
minutes after 16:00 UT 



90 




20 30 40 50 60 70 80 90 
minutes after 16:00 UT 



70 80 90 
:00UT 



Fig. 4. — Evolution of coronal-averaged temperature (a), density (b), and pressure (c) of a flux tube 
with varying fraction of energy flux carried by beam driven upflows, T/H = 0-0.8. Also plotted in 
dashed lines are arbitrarily scaled total heating rate in panel (b) and observed UV light curve in 
panel (c). 
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Fig. 5. — The differential emission measure (DEM) derived from multiple flare loops from 16:41 to 
16:43 UT, with different 7^ but the same total heating rates. 
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Fig. 6. — Left: the sum of the total heating rate (purple), ad hoc heating rate (red), and beam 
heating rate (blue) in all flare loops, in comparison with the observed total UV light curve (black) 
and HXR 25-50 keV light curve (orange), both arbitrarily scaled. Right: peak temperature and 
density distributions in over 5000 modeled flare loops with lengths in the range of 35-55 Mm and 
cross-sectional area of 1" by 1". 
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Fig. 7. — Time evolution of the coronal DEM from the model of the 2005 May 13 flare. The 
high emission measure at low temperature at the beginning of the flare is produced by background 
heating. 
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Fig. 8. — Left: Comparison of the computed (red) and GOES observed (black) SXR fluxes at 
0.5-4 A (dotted) and 1-8 A (sohd). 
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Fig. 9. — Comparison of the synthetic SXR hght curves from the model (red) with RHESSI ob- 
servations (black) in the 3-6 keV, 6-12 keV, and 12-25 keV bands. The vertical dotted lines and 
letters a-e indicate the times when the spectral comparison is plotted in Figure [TOl The discon- 
tinuities in the light curves at 47 min and 57 min are due to changes of the RHESSI attenuator 
state, when the thick shutter is inserted and later removed. 
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Fig. 10. — Comparison of synthetic SXR spectra from the model (red) with the RHESSI ohservation 
(black) at a few times during the flare. The time of each panel is indicated by vertical dotted lines 
in Figure m For a complete view, a movie is available lonlinel which compares the observed (black) 
and computed (red) spectra from 16:37 to 17:15 UT, with 20 s cadence. 
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Fig. 11. — Top: time profiles of the ratio (black) of the modeled SXR spectrum to that observed 
by RHESSI averaged over the 6-15 keV band, and slope (red) of the ratio versus photon energy in 
the 6-15 keV band. Bottom: time histograms of the number of flaring pixels identified from the 
UV foot-point emission. 



36 



Q 



5x10" 
4x10' 
3x10" 



M 2x10' 



> 

3 



1x10" 



obs. UV flux 
comp. UV flux 




20 30 40 50 60 70 80 90 
minutes after 16;00 UT 



6.0 














^M\ obs. total flux 


^ 2.5 


- 






^ 2.0 


_ 


• 

1 


\ comp. decay flux \ 


o 




1 


\ 


zz, 1.5 


- 


/ 


V 


X 




1 


N 


^ 1.0 


- 


4 
1 


- N 


^ 0.5 


- 


1 


F ^^^^^^^^^^^ 


0.0 




• 


y. T~"~~~~" 



20 30 40 50 60 70 80 90 
minutes after 16:00 UT 



Fig. 12. — Comparison of the synthetic transition region C iv hne emissions (red) in the decay 
phase with the UV flux observed in the TRACE 1600 A band (black) for one flare loop (left) and 
for all flare loops (right). In the right panel, the observed decay flux is derived by subtracting from 
the observed light curve the full Gaussian profile that represents the UV rise. 



